Multicomponent Skyrmion lattices and their excitations 
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We study quantum Hall ferromagnets with a finite density topologically charged spin textures in the presence 
of internal degrees of freedom such as spin, valley, or layer indices, so that the system is parametrised by a 
d-component complex spinor field. In the absence of anisotropies, we find formation of a hexagonal Skyrmion 
lattice which completely breaks the underlying SU (d) symmetry. The ground state charge density modulation, 
which inevitably exists in these lattices, vanishes exponentially in d. We compute analytically the complete low- 
lying excitation spectrum, which separates into d 2 — 1 gapless acoustic magnetic modes and a magnetophonon. 
We discuss the role of effective mass anisotropy for S(7(3)-valley Skyrmions relevant for experiments with 
AlAs quantum wells. Here, we find a transition, which breaks a six-fold rotational symmetry of a triangular 
lattice, followed by a formation of a square lattice at large values of anisotropy strength. 
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Introduction. Skyrmions in quantum Hall ferromagnets fljj] 
present an early example of topological excitations in con- 
densed matter physics, a topic of persistently large interest, 
which by now has grown to encompass itinerant magnets with 
spin-orbit coupling |2j,[3[] and spinor condensates [4]. In the 
quantum Hall effect (QHE), these topological textures carry a 
quantised charge, endowing them with stability and providing 
possible probes via Coulomb interactions and charge trans- 
port. The low-energy modes of these defects have possibly 
been observed in NMR measurements [5]. The textures are 
also exceptionally tunable, as their density, and hence the rel- 
ative strength of interactions, can be modified by a gate po- 
tential or magnetic field. The possibility of controlling spin 
by coupling to electron charge makes these systems interest- 
ing candidates for spintronics applications. 

In this work, we discuss the physics of Skyrmions with an 
enlarged internal space. This space minimally consists of two 
components which are, most simply, represented by physical 
up and down spins - but other degrees of freedom frequently 
play a role. Early examples are the layer index in double quan- 
tum wells of GaAs, where interlayer phase coherence sponta- 
neously develops , or the valley degeneracy in semicon- 
ductors [8]. Recently, there has been a further proliferation, in 
the form of spin, layer and valley degeneracies in graphene 
ii, and cold atoms 0], where SU(N) internal degrees of 
freedom and Skyrmions have emerged as interesting topics. 

At finite density of defects and in presence of Coulomb in- 
teractions, Skyrmions can form a crystal (SC), analogous to 
Skyrme crystals of nuclear physics fioll . The properties of 
SCs have been extensively discussed in the SU (2) case, and a 
large number of results in a wide variety of settings have been 
obtained using Hartree-Fock (HF) theory for SU (4) systems 
111 111 . Most of these studies are based on the numerical solu- 
tion of the HF equations of motion, being hard to perform due 
to presence of several lengthscales. 

This motivates us to develop an analytical approach to 
the general case of Skyrmion crystals with an internal space 
parametrised locally by a d-component complex spinor. The 



simplifying feature of our model is that the spinors live in 
the complex projective space CP d_1 , which supports exact 
multi-Skyrmion solutions irrespective of the value of d, and, 
most crucially, the assumption of large value of spin stiffness. 

Our programme consists of two parts. First, we derive a set 
of largely analytical results on the case of Skyrmion crystals 
which are described locally by SU (d) spinors. We find that 
macroscopic degeneracy, present in the noninteracting case 
JI2I1 . is lifted in favour of a hexagonal SC with the amplitude 
of charge density modulations decaying exponentially with d. 
The excitation spectrum, obtained by studying small pertur- 
bation of the crystal, exhausts the number of gapless modes 
consistent with the original degeneracy. Here, d 2 — 1 of these 
modes correspond to generators of the SU (d) Lie algebra, in 
accordance with the Goldstone theorem, supplemented by a 
magnetophonon excitation, appearing as a result of broken 
translational invariance. The magnetophonon has a dispersion 
reminiscent of a 2D Wigner crystal in a magnetic field [13]. 

Second, as an application of our theory, we consider explic- 
itly a model of multi-component quantum Hall ferromagnet, 
motivated by experiments on semiconductor quantum wells. 
The effects of anisotropies, present in these systems, which 
break SU(d) symmetry, are studied. With increasing d, the 
number of possible types of anisotropies grows quickly along 
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FIG. 1: (Color online). Positions of zeros of the basis theta-functions 
8 p (z) within a single unit cell of a Skyrmion lattice in the SU(3) 
case. The lattice vectors are indicated by arrows. 
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with the d 2 — 1 generators of SU(d), and we concentrate on 
the simplest nontrivial, yet instructive, and hitherto unstudied 
case of SU(3) Skyrmions. These are of experimental inter- 
est for quantum wells grown in a [111] direction, in particu- 
lar a nascent effort on AlAs, where one finds the appearance 
of three almost degenerate valleys [8]. The latter are subject 
to a 'nematic' anisotropy recently discussed in the context of 
SU(2) valley degrees of freedom 11411 . which we generalise to 
the case of an arbitrary number of valleys. We show how this 
anisotropy changes the nature of the crystalline ground state 
by combining spontaneous symmetry breaking with the exter- 
nally imposed one, leading to a collection of different crystal 
structures. We also determine the gaps appearing in the low- 
energy excitation spectrum. 

The remainder of the paper is organised as follows. First, 
we introduce the CP d ~ 1 model describing spatial variations 
of internal degrees of freedom, taking into account interac- 
tions between electrons. An exact HF ground state in terms 
of spinors living in an analytic subspace is presented. Start- 
ing from this state we derive linearized equations of motion 
describing excitations of Skyrmion lattices. A generalization 
of Bloch theorem leads, in the isotropic case, to a reduction 
of a system of 2d 2 equations into d 2 decoupled pairs of linear 
equations, which analytical solution provides the spectrum to- 
gether with the eigenfunctions of the excitations. Next, we 
consider a model of semiconductor, possessing three degener- 
ate valleys, which takes into account the 'nematic' anisotropy, 
and present results for the phase diagram and excitations. 

The model. We consider a quantum Hall ferromagnet at fill- 
ing factor v ~ 1 with internal degrees of freedom described by 
a d-component spinor w n (r), which corresponds to the Slater 
determinant state \^>) = II r [J^ ^ iy n (r)cjj(r)]|0), where the 
spinor is normalized as $D_ =1 |^n(r)| 2 = 1 at every point r, 
& n (r) creates an electron in the n-th spin component in the 
lowest Landau level (LLL) eigenstate with position r. The 
latter are given in the coherent state representation, and |0) is 
the electron vacuum. For example, in the case of a bilayer 
quantum Hall system, components of the spinor correspond 
to layer and spin indices, and in graphene to spin and valley. 
The Hamiltonian of the system is given by two contributions 
H = Hcp + Hi n t, where Hep is proportional to the spin stiff- 
ness ps, and H- m t describes Coulomb interactions. A generic 
form of Hep, invariant under SU(d) rotation and containing 
only first order gradients reads J3.Q 



in the lowest order of the gradient expansion. Here Q(r) is 
the topological charge density, which is proportional to the 
electronic charge density in the quantum Hall effect 111 Il5ll . 
and is related to a Pontryagin index of the normalized spinor 
field ill 211 . In the following we consider textures in the form 
of Skyrmion crystals and fix the topological charge, given by 
the integral of Q(r) inside the unit cell (u.c), to 2nd. The 
two lattice vectors of the unit cell, which define a parallelo- 
gram with opening angle 0, can be chosen as 71 = %\[d and 
72 = irVdr, where r = e 1 ^ (real/imaginary parts of 712 cor- 
respond to x, y components of the vectors). We measure the 
lengths in units of 21b / — l|ImT. 

The energy Hqp is minimized by spinors, with entries ana- 
lytic functions of z = x+iy, whose number of zeros inside the 
u.c. fixes the topological charge sector. Without interactions 
there is a macroscopic degeneracy coming from the arbitrari- 
ness in the choice of positions for these zeros. In presence 
of £/i n t this degeneracy is lifted, which can lead to formation 
of a Skyrmion crystal. It is impossible to write an analytic 
expression for the ground state in the interacting case, thus 
we propose a variational approach in which the trial state is 
still expressed in terms of analytic functions w a (z). Although 
the Euler-Lagrange equations, derived from Hep, are highly 
nonlinear, our procedure is reminiscent of the LLL projection, 
familiar from the theory of Abrikosov vortex lattices and the 
QHE physics. Our conjecture is that the zero modes of the 
Hessian of Hqp, spanning the analytic subspace, are sepa- 
rated from the higher (non-analytic) ones by a spectral gap of 
the order p$, because the average topological charge density, 
analogous to a magnetic field strength, is independent of d. 
Further, the energy scale associated with the Coulomb inter- 
action is psy/nslB, where rig is the Skyrmion density and Is 
describes the magnetic length. This scale is smaller than p s 
provided that ngl 2 B <C 1, i.e. if the filling factor remains close 
to one. This justifies our variational treatment. 

The method. We choose a basis of analytic functions de- 



fined by their quasi-periodicity 9(z +7) = e a i z + b -r 
rii 71 + ri272, where a 7 , b 7 are complex numbers and m,2 € 
Z. Using the properties of these functions under transla- 
tions f 12/d 6 p (z) = e-^ d+i y% +1 {z) and f ll/d 6 v (z) = 
e 2 ™P/ d O p {z), we obtain the basis of linearly-independent 
theta-functions with a given topological charge inside the u.c. 



= i7T-rd(n-§)(ri-l-§)+2i\/d(n-§)z 



(3) 



H C p = 2 Ps / d 2 r 



.(Vw,Vw) (w, "Vw)(Vw, w) . 



(1) 



(w,w) (u>,u>) 2 
where the integral is taken over the two dimensional plane, 
and the scalar product is defined as (a, b) = J2 n a *i( r )^«( r )- 
The Hamiltonian has a local "gauge" symmetry, in other 
words, multiplying a spinor by an arbitrary function, w a (r) = 
h(r)w a (r), does not affect the energy. The electron interac- 
tions with a potential V(r — r') are described by 



Hurt = \ J d 2 r d 2 v' V(t - r')Q(r)Q(r'), 



(2) 



where p runs from to d — 1, in agreement with Riemann- 
Roch theorem lllal . From the numerical minimization proce- 
dure, which was further supported by a linear stability anal- 
ysis, we find that the energy minimum is achieved for a 
spinor formed from basis functions with equal amplitudes (up 
to a global SU(d) rotation) w a (z) = (9q, . . . , 8d-i), and 
r = e M / 3 . This corresponds to a charge distribution with the 
full symmetries of a hexagonal lattice, although the individual 
spinor components possess a lower symmetry. 

In order to calculate the excitations spectra of the SC we use 
the basis formed from a set of quasi-periodic analytic func- 
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tions Xpk( z ) obtained from 9 p (z) by translations by vector k 
with components 2-Km\jN\ and l-Km^l^i, where mi, £ 
Z and Ni and N2 is the number of unit cells on a torus. An 
arbitrary analytic spinor can be written in terms of Xpk( z ) as 

^(^Ewp^EL^W^W' o = 0..d-l (4) 



see 1 17] and a related work JT§], where an analogous approach 
was formulated to a (7(1) case for the problem of vortex lat- 
tices in Bose-Einstein condensates. Here we introduce arbi- 
trary complex square matrices U of coefficients with linear 
dimension d. We consider time-dependent HF, with dynam- 
ics generated by a Lagrangian given as a sum of the energetic 
term H lnt and a Berry's phase term, which is first order in 
time derivatives of the spinor components. Variation of the 
Lagrangian with respect to parameters U a b (k) gives a system 
of coupled equations describing the motion of a Skyrmion lat- 
tice, constrained to the analytic subspace 



d 2 r Xbk( z ) r-d^q , (w, d t w) 



Wa(z) 



(w, w) dt (w, w) 

-4w a (z) J A T V(r - r')Q(r')d 2 r'} = 0. (5) 

Equation (O represents one of the main results of our paper. 

LSWT and Block theorem. The system of equations (0 is 
very general and is not amenable to analytic treatment, mostly 
because the velocity field dw a /dt is not explicit. However, 
due to the high degree of symmetry of our reference SC state, 
it is possible to obtain full analytical results at arbitrary d for 
small perturbations of the lattice, i.e. for linearized equations 
of motion. The matrices of coefficients now can be written as 
a sum W a b(k) = <5 a fc<5k.o + "afc(k), where w a &(k) are small. 
The quasi-periodicity of theta-functions and the periodicity 
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FIG. 2: Excitation spectra of a SU (3) Skyrmion lattice, in the case 
of Coulomb interactions, along a diagonal direction in the Brillouin 
zone, fci = fe. The lowest lying mode at small fei is a mag- 
netophonon. All remaining modes are linear gapless acoustic ex- 
citations, (inset) appearance of gaps in the SU(3) case at finite 
anisotropy strength. 



(with elementary periods ji/d and 72 / d) of the topological 
charge density of the SC leads to a generalized Bloch theo- 
rem, and the system of equations fully decouples into d 2 pairs 
of linear equations. These are readily solved, resulting e.g. in 
the excitation spectrum for a SU(3) case presented in Fig. 2. 

Results. The magnetophonon mode at small energies has a 
dispersion |k| Q , with a ~ 3/2 in the case of Coulomb inter- 
actions, like the dispersion of a 2D Wigner crystal in a mag- 
netic field. For a short-range interaction a ~ 2, while the 
spectrum is linear at small energies, a ~ 1, for logarithmic 
(2D Coulomb) interactions. In the large-c? limit, topological 
charge density variations of a SC have an exponentially small 
amplitude, for example for a square lattice we find 

Q(r) = 2/Tr-Me-^[cos(2Vdx)+cos(2Vdy)]+0(e-* d ), 

Coulomb interactions suppress charge density fluctuations, 
and, for a large number of components, it becomes possible to 
arrange for a nearly-uniform twist of the spinor. Interestingly, 
similar questions arise in the studies of SU(N) superconduc- 
tors J3. 

SU(3) case. Let us discuss application of our general meth- 
ods to a potentially realizable case of Skyrmion crystals 
in AlAs quantum wells, where the spin degree of freedom 
is quenched due to applied magnetic field, and the valley- 
degeneracy leads to the three-component spinors with coeffi- 
cients Oj, which we describe in zeroth order by a SU (3) sym- 
metric theory. Effective mass anisotropy within the valleys 
with relative angle ip generates a symmetry-breaking term in 
the Hamiltonian describing a 'nematic' anisotropy 



-A c 1 +2A k 



Mi 



<d 2 |a,- 



(6) 



where Aq 



ie 2 /sIb is the exchange energy in the 



isotropic case, and Ib = WHc/eB is the magnetic length. 
With K(x) denoting the complete elliptic integral of the first 
kind, and n = (ci — C2)/2, the coefficients are given by 



2 K(l- 1/A 2 ) 
ci = 7= , c 2 



2 K{1 - b 2 /a 2 ) 



(7) 



where a 2 — A cos 2 ^ip + A 1 sin 2 ^<p and b 2 = A sin 2 ^ip + 
A -1 cos 2 \<p>, and A = {nix/my) 1 / 2 stands for the anisotropy 
parameter. Notice that x > 0, describing repulsive interac- 
tion between the different components a, of the spinor. This 
reduces to the form derived in IU4I1 for the SU(2) case. In 
the isotropic limit, rn x = m y , the ground state, for a given 
filling factor, is a triangular Skyrmion lattice, which trans- 
forms into a square one with increasing value of the anisotropy 
strength, see Fig. 3. This transition is due to the relative scal- 
ing of Coulomb and nematic energies with Imr. We do not 
find the nematic phase transition of [14] even at large values 
of anisotropy strength. Note that although valley polarization 
is favoured by the anisotropic term, its Coulomb energy cost 
turns out to be so large that the different spinor amplitudes 
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FIG. 3: a) Angle (j> in the 5(7(3) case as a function of nematic 
anisotropy strength. The crystal is hexagonal in the isotropic case 
and becomes square at large Jjv = xexp(7rd/2)/7r\/27r|^ — 1| 3 ^ 2 . 
b) Scaling of gaps with nematic anisotropy strength. The first gap 
shows linear behavior, while the other two increase as a square root, 
(upper inset of Fig.a) Schematic picture of the valleys in AlAs quan- 
tum wells in the [111] direction showing the three-fold valley de- 
generacy. The effective mass within each valley is anisotropic with 
m x /m y w 4. (lower inset of Fig.a) topological charge densities of a 
Skyrmion crystal from weak to strong anisotropy. 



in our variational ground-state remain equal, independently of 
the value of x. 

Finite gaps appear in the spectrum in presence of 
anisotropy. The first three modes remain gapless; they corre- 
spond to magnetophonon and to the two remaining continuous 
symmetry generators associated with a two-dimensional Car- 
tan subgroup of SU(3). The other six modes are now gapped, 
so that the collective spectrum forms the pattern shown in 
Fig. 3. The first two gaps scale linearly with the anisotropy 
strengths, while the behavior is ~ \/Jn for the other four. 

Energy scales. In the SU(3) valley-degenerate case we ob- 
tain at small anisotropies An = Aox/27rZf; ~ jfg (A — 

l) 2 , and p s = e 2 /16V^el B ~ 5.2 K, with A - 131 K. 
The anisotropy energy Aqx ~ 2.5 K is much smaller than 
the spin stiffness scale Airps ~ 65 K, where we assumed e = 
10eo for AlAs at v = 1 and electron density 2.5 x 10 11 cm _1 



j20ll . The characteristic Coulomb energy is ~ ps\v — ll 1 / 2 
and can be made small compared to the stiffness by changing 
the filling factor. For relatively small values of \v — 1| ~ 0.1 
we estimate the value of dimensionless parameter Jn ~ 13.5 
which corresponds to a square lattice, see Fig. 3. 

Conclusions. We have studied SU(d) quantum Hall 
ferromagnets with finite density of topological defects 
(Skyrmions). In presence of Coulomb interactions the ground 
state of the system is a triangular Skyrmion lattice, which 
breaks all the symmetries of the internal SU(d), as well 
as translational symmetry, generating d 2 Goldstone modes, 
d 2 — 1 of which correspond to breaking of the SU(d), and 
the remaining one being a magnetophonon mode. We have 
explored the phase diagram of a SU(3) QHE ferromagnet, 
which is expected to be relevant to semiconductor nanostruc- 
tures with valley-degeneracies [8]. The lattice tilts continu- 
ously as a function of anisotropy strength and becomes square 
through a phase transition at a critical anisotropy value. In 
this case, three of the excitation branches remain gapless, 
while others acquire gaps exhibiting different scaling with 
anisotropy strength. 

It would be interesting to extend our theory to the case 
of large, but finite, spin-stiffness, and to derive an effective 
sigma-model for the Skyrmion lattice in the projected sub- 
space. The difficulty which one immediately faces is that 
long-wavelength SU(N) rotations are not preserved under 
projection onto the analytic subspace. Another important 
question is related to entanglement of internal degrees of free- 
dom, as studied for simple Skyrmions structures in 112 ill . 
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